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ABSTRACT 

Supermassive black holes (SMBHs) have been detected in the centers of most nearby massive galax- 
ies. Galaxies today are the products of billions of years of galaxy mergers, but also billions of years of 
SMBH activity as active galactic nuclei (AGNs) that is connected to galaxy mergers. In this context, 
detection of AGN pairs should be relatively common. Observationally, however, dual AGN are scant, 
being just a few percent of all AGN. In this Letter we investigate the triggering of AGN activity in 
merging galaxies via a suite of high resolution hydrodynamical simulations. We follow the dynamics 
and accretion onto the SMBHs as they move from separations of tens of kiloparsecs to tens of parsecs. 
Our resolution, cooling and star formation implementation produce an inhomogeneous, multi-phase 
interstellar medium, allowing us to accurately trace star formation and accretion onto the SMBHs. We 
study the impact of gas content, morphology, and mass ratio, focusing on AGN activity and dynamics 
across a wide range of relevant conditions. We test when the two AGN are simultaneously detectable, 
for how long and at which separations. We find that strong dual AGN activity occurs during the late 
phases of the mergers, at small separations (<1-10 kpc) below the resolution limit of most surveys. 
Much of the SMBH accretion is not simultaneous, limiting the dual AGN fraction detectable through 
imaging and spectroscopy to a few percent, in agreement with observational samples. 
Subject headings: galaxies: nuclei — galaxies: interactions — galaxies: active 



1. INTRODUCTION 

The exceptional spatial resolution of the Hubble Space 
Telescope and of Chandra led to the exciting discovery of 
the first three cases of dual active galactic nuclei (AGNs) 
in the center of the same galaxy, in LBQS 0103-2753, 
NGC 6240 and Arp 299 (Junkkarinen et al. 2001; Ko- 
mossa et al. 2003; Ballo et al. 2004). Recently, more ob- 
servations were focused on detecting spatially resolved 
dual AGN via various techniques (Gerke et al. 2007; 
Smith et al. 2010; Hudson et al. 2006; McGurk et al. 
2011; Rodriguez et al. 2006; Bianchi et al. 2008; Barth 
et al. 2008; Comerford et al. 2009b; Liu et al. 2010; Pi- 
concelli et al. 2010; Comerford et al. 2009a; Green et al. 
2010; Liu et al. 2011; Fu et al. 2011; Koss et al. 2010, 
2011; Fabbiano et al. 2011). 

If most galaxies host a SMBH (Ferrarese & Ford 2005) 
and galaxy mergers trigger quasar activity (e.g., Di Mat- 
teo et al. 2005) then one expects that dual AGNs should 
be common. Observationally, however, AGN pairs are 
rare (at most a few percent, see references above). The 
prediction of the timescale on which SMBH pairs can 
be observed as double quasars (Hennawi et al. 2006; Liu 
et al. 2011, and references therein) is a key diagnostic of 
SMBH merger rates (of paramount importance for grav- 
itational wave searches) and AGN triggering. 

Foreman et al. (2009) discuss how the relationship be- 
tween the lifetime of an active SMBH, £agn, and the 
merging timescale, £ m erg, plays a fundamental role in de- 
termining the observability of AGN pairs. If one assumes 
that most galaxies host SMBHs, that AGN/quasar activ- 
ity is triggered by galaxy mergers, and that the lifetime of 
quasars equals the merger timescale, £agn — ^merg, the 
probability of observing a dual quasar should be close 



to unity, if we do not consider additional factors, such 
as obscuration. If £agn ^ ^merg (Di Matteo et al. 2005; 
Hopkins et al. 2005), and/or if there is a delay in the trig- 
gering of the two quasars, then one might have ceased its 
activity before the other started. Foreman et al. (2009) 
also notice that the distribution of physical separations 
for luminous quasar pairs in the Sloan Digital Sky Survey 
peaks below 30 kpc, the lower limit of the physical reso- 
lution that can be resolved in the survey. The paucity of 
optically selected quasar pairs on galactic scales (~ 0.1% 
at L > 10 45 erg s _1 , Foreman et al. 2009; Hennawi et al. 
2006) points toward non-simultaneous activity at large 
separations. 

At lower levels of activity, Comerford et al. (2009b); 
Liu et al. (2010); Shen et al. (2011) find that about 2-5% 
of optically-selected (via K/3 and [OIII] lines) AGNs are 
in pairs, and 30% off-set from the host center (hinting 
at inspiral). On the other hand, about 30% of X-ray de- 
tected sub-mm galaxies at z « 2 are in pairs (Alexander 
et al. 2003), suggesting that the occurrence of AGN pairs 
is possibly underestimated by optical selection. However, 
some observed AGN pairs are likely to be caused by gas 
kinematics rather than true dual AGN (Shen et al. 2011; 
Fu et al. 2012). 

Volonteri et al. (2009) and Yu et al. (2011) also discuss 
models that reconcile theoretical merger rates of SMBHs 
and galaxies with the small fraction of binary quasars 
and dual AGN respectively. The lifetime of AGN, the 
gas content of the host galaxies and the dynamics of the 
merger are the main factors that explain the paucity of 
observed AGN pairs. Understanding the occurrence of 
AGN pairs therefore requires a thorough investigation 
and understanding of the detailed physical conditions de- 
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scribing the evolution of SMBHs during mergers. 

In this Letter we investigate the theoretical expecta- 
tions for detections of dual AGNs. We use high reso- 
lution simulations to study the dynamics, statistics, and 
observability of AGN pairs to test if the scarce number of 
quasar /AGN pairs is consistent with limited observabil- 
ity of simultaneous activity of SMBHs involved in galaxy 
mergers. 

2. SIMULATIONS 

We investigate the dynamical and accretion history 
of SMBHs in merging galaxies via a suite of very high 
resolution (< 20 pc) smoothed particle hydrodynamics 
(SPH) simulations. We present here two new simulations 
with mass ratios 1:2, one between two gas-rich spiral 
galaxies and the other between a gas-poor elliptical and 
a gas-rich spiral. We supplement our analysis with a 1:10 
mass ratio simulation between gas-rich spirals from Cal- 
legari et al. (2011). This suite focuses on galaxy mergers 
most relevant to the build up of a galaxy in the ACDM 
cosmology. As mergers between equal mass galaxies are 
rare, we probe the range of mergers that are more com- 
mon but are expected to lead to efficient SMBH pairing. 
The peak of the cosmic merger rate occurs at higher red- 
shift, so we begin our simulations at z = 3. We include 
the elliptical-spiral merger to study the impact of gas 
fraction and morphology on SMBH pairing and accre- 
tion. This range of parameters allows us to study SMBH 
and host dynamics and co-evolution across a range of 
cosmologically relevant merger conditions. 

The methodology of how our spiral galaxies are initial- 
ized, as well as the details of the 1:10 merger, is described 
in Callegari et al. (2011). We briefly summarize here the 
galaxy models that we use in the 1:2 mergers. The spiral 
galaxies contain a dark matter halo, a disk composed of 
stars and gas, and a stellar bulge. The halo is represented 
by a spherical Navarro et al. (1996) profile with spin pa- 
rameter A = 0.04. An exponential disk of stars and gas 
is included with total mass 0.04M v i r and gas fraction 
f g = 0.3. We note that observations of more actively star 
forming galaxies at high redshift suggest even higher gas 
fractions (Tacconi et al. 2010). The stellar bulge is repre- 
sented by a spherical Hernquist (1990) model with mass 
0.008M v i r . The primary galaxy in the 1:2 and 1:10 spiral- 
spiral mergers and the secondary galaxy in the elliptical- 
spiral merger have M vir = 2.24 x lO n M . The elliptical 
galaxy consists of a dark matter halo and a stellar com- 
ponent, each represented by a Hernquist (1990) profile. 
The halo has spin parameter A = 0.04 and scale length 
22 kpc, chosen to resemble a Navarro et al. (1996) profile 
following the method in Springel et al. (2005), assuming 
concentration c = 3. The stellar component has total 
mass 0.05M v i r , no rotation, and scale length 0.5 kpc, 
chosen based on the scale length of early type galaxies 
from SDSS (Shen et al. 2003). 

In both 1:2 mergers represented here, dark matter par- 
ticles have masses of 1.1 x 1O 5 M0 and softening length 
30 pc, gas particles have masses of 4.6 x 1O 3 M and 
softening length 20 pc, and star particles have masses of 
3.3 x 1O 3 M and softening length 10 pc. This resolution 
allows us to consistently track the dynamics and evolu- 
tion of the SMBHs down to the formation of a SMBH 
pair while minimizing their excursion from the centers 
of their galaxies. The high resolution also enables us to 



track the dynamics of gas in the galaxies, leading to accu- 
rate determinations of SMBH accretion rates, feedback, 
and star formation during different stages of the merger. 
Black hole masses are chosen to be consistent with the 
M B n -M bu i ge relation (Marconi & Hunt 2003). A single 
particle representing the SMBH is placed at the center of 
each galaxy. The initial SMBH masses are 4 x 1O 7 M and 
3 x 1O 6 M in the elliptical-spiral merger and 3 x 1O 6 M 
and 1.5 x 1O 6 M in the 1:2 spiral-spiral merger for the 
primary and secondary SMBHs, respectively. For refer- 
ence, the SMBH masses in the 1:10 merger from Callegari 
et al. (2011) are 6 x 1O 5 M and 6 x 1O 4 M . 

Each simulation begins with the galaxies at a separa- 
tion equal to the sum of their virial radii. Orbital pa- 
rameters were chosen based on results from cosmological 
simulations (Benson 2005). In the 1:2 elliptical-spiral 
merger, the galaxies have initial eccentricity 0.98 and 
pericentric distance 19% of the virial radius of the pri- 
mary galaxy. In the 1:2 spiral-spiral merger, the galaxies 
have initial eccentricity 1.02 and pericentric distance 26% 
of the virial radius of the primary galaxy. All mergers 
are planar and prograde. We note that this orientation 
may produce stronger gas inflow than inclined mergers, 
likely maximizing SMBH accretion and, therefore, dual 
SMBH activity. 

We performed our simulations using GASOLINE, a 
SPH AT-body Tree code (Wadsley et al. 2004; Stadel 
2001). GASOLINE includes a physically motivated 
prescription for star formation and supernova feedback 
(Stinson et al. 2006), as well as a recipe for black hole 
physics (Bellovary et al. 2010). We have chosen our pa- 
rameters 1 to match those in the literature for simulations 
of comparable resolution , which are able to produce re- 
alistic galaxies (e.g. Governato et al. 2010). Stars are 
able to form if the parent gas particle reaches a threshold 
density of 100 amu cm -3 and is below a temperature of 
6000 K. Supernovae release Es n — 10 51 erg into the sur- 
rounding gas based on the blastwave formalism of Stinson 
et al. (2006). In order to prevent an unphysical burst of 
supernovae as the simulation begins, we relax the galax- 
ies in isolation, gradually increasing the star formation 
efficiency, c*, from 0.005 to 0.015 over 10 8 years. Black 
holes accrete gas through Bondi-Hoyle accretion. This 
accretion gives rise to thermal feedback, which we model 
as E = 6f€ r Mc 2 (Springel et al. 2005), with radiative ef- 
ficiency e r = 0.1 and feedback efficiency e/ = 0.001 as in 
Bellovary et al. (2010). The feedback energy is imparted 
on the nearest gas particle to the black hole. 

3. RESULTS 

In Table 1, we present the results of our simulations. 
We determine the fraction of the total simulated time 
that each SMBH is active, given a threshold for observ- 
ability. We also show the time spent as a dual AGN, 
when both SMBHs exceed the given threshold at the 
same time, expressed as a fraction of the total time that 
one or both SMBHs are active at that threshold. The 
dual AGN fraction is also given with additional con- 
straints to mimic observational limitations. The sep- 
aration thresholds represent possible spatial resolution 

1 The new simulations presented in this Letter use an updated 
version of GASOLINE. We refer the reader to Callegari et al. (2011) 
for a list of parameters used in the 1:10 simulation. 
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Table 1 

BHi and BH2 represent the percentage of the total simulated time that the primary and secondary SMBHs are active, respectively, above 
the given activity threshold (L42 = 10 42 erg s _1 ; L43 = 10 43 erg s _1 ; L44 = 10 44 erg s _1 ; accreted mass at each timestep is converted 
into a bolometric luminosity assuming a radiative efficiency of 10%). The total simulated time for each merger is: 1.35 Gyr for the 1:2 
Spiral-Spiral simulation, 1.1 Gyr for the 1:2 Elliptical-Spiral simulation, and 2.56 Gyr for the 1:10 Spiral-Spiral simulation. The Dual 
AGN column gives the percentage of time that both SMBHs are active (over the total time that one or more SMBHs are active at the 

given threshold). The remaining columns give the dual activity percentage with additional constraints, based on imaging or spectroscopic 
detect ability, d represents the absolute separation between the SMBHs. Av is the absolute velocity difference between the SMBHs. 




t (Gyr) 



Figure 1. Evolution of the 1:2 spiral-spiral merger. Upper panel: 
separation of the two black holes as a function of time. The inset 
shows a zoomed-in view of the separation at late times. Middle 
panel: star formation rate from both galaxies as a function of time. 
Lower panel: bolometric luminosities of the two black holes as a 
function of time. The blue, dotted line represents the primary 
SMBH. The red, solid line represents the secondary SMBH. 

cutoffs for separating the AGN with imaging. Spectro- 
scopic duals require sufficiently large velocity offsets that 
two sets of emission lines are discernible. We note that 
these cutoffs are absolute, not projected, quantities in 
our analysis; they are meant to show the qualitative ef- 
fects of observational limits. 



1:2 Spiral-Spiral merger. We show in Fig. 1 the evo- 
lution of the SMBH separation, host galaxy star forma- 
tion, and SMBH accretion. Star formation rates include 
stars formed in both the primary and secondary galax- 
ies. Before the second pericenter passage, the SMBHs 
accrete relatively little, increasing in mass by 20%-30% 
over a Gyr of evolution. Accretion is not well correlated 
between the SMBHs, and bolometric luminosities gener- 
ally remain at < 10 43 erg s _1 . Following the second peri- 
center passage, however, tidal torques concentrate gas in 
the centers of the two galaxies, leading to efficient SMBH 
accretion accompanied by strong star formation, which 
may complicate detection of the AGN. Over the next 200 
Myr, the SMBHs increase in mass by factors of 30%-60%. 
During this phase, when the SMBHs are separated by 
less than 10 kpc, the AGN bolometric luminosities reach 
10 44 erg s _1 , and the activity is better correlated, trig- 
gering after each pericenter passage. These luminosities 
are sustained for ~ 35 Myr by each SMBH separately, 
and for ~ 10 Myr simultaneously. We end this simula- 
tion at t = 1.35 Gyr, as the SMBHs have reached our 
resolution limit at separations of 10 pc. 

Fig. 2 2 shows the amount of time both SMBHs are 
active at a given separation and luminosity ratio for dif- 
ferent thresholds. At low luminosity, the SMBHs are 
active for 60%-70% of the simulation and the observ- 
ability of dual AGN activity traces the orbit of the host 
galaxies. There is significant activity at large separa- 
tions as the galaxies spend most of their time at or near 
apocenter. However, a higher activity threshold selects 
for the phases of the merger where AGN triggering is 
strongest. The dual AGN activity at higher thresholds 
occurs mainly at small separations (< 10 kpc), follow- 

2 Additional figures for the other simulations 
and with log-scale separations are available at 
http : / / www. astro . lsa. umich . edu / ~ s vanwas / dual AGN . html 
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Figure 2. Dual AGN observability timescale at a given SMBH separation, activity threshold, and luminosity ratio or velocity offset for 
the 1:2 spiral-spiral merger. In each case, the threshold is imposed upon both SMBHs. 



ing the second and subsequent pericenter passages. We 
show also the distribution of velocity offsets between the 
SMBHs during dual AGN activity. The longest episode 
of dual activity above 10 44 erg s _1 occurs following the 
third pericenter passage, below separations of ~ 2 kpc 
(see inset of Fig. 1), near apocenter. The velocity offset 
is small, therefore, and we find that c± 75 percent of the 
dual activity is at Av < 150 km s _1 . In 1:1 mergers, Di 
Matteo et al. (2005) and Hopkins et al. (2005) also find 
that the strongest SMBH activity occurs at small sepa- 
rations, although they do not distinguish between single 
and dual AGN activity. 

1:2 Elliptical-Spiral merger. In the 1:2 elliptical- 
spiral merger, the elliptical galaxy initially contains no 
gas. However, gas stripped from the companion galaxy 
during the first and second pericenter passages cools and 
settles at the center of the elliptical, allowing the pri- 
mary SMBH to begin accreting at around t = 0.7 Gyr 
(Fig. 3). Since the primary SMBH is an order of mag- 
nitude more massive than the secondary, it produces a 
higher luminosity once gas is available. As in the spiral- 
spiral merger, the strongest SMBH accretion occurs fol- 
lowing the second and subsequent pericentric passages. 
Above 10 44 erg s _1 dual AGN activity occurs for ~ 9 
Myr, mostly under 10 kpc separations. We stop this 
simulation at t = 1.1 Gyr. At this time, the stars and 
gas surrounding the secondary SMBH have been tidally 
disrupted, preventing further AGN activity. The SMBHs 
remain at a separation of hundreds of pc. The dynamical 
friction timescale at this separation, ignoring the effects 
of gas, is of order a few hundred Myr (Colpi et al. 1999), 
therefore the two SMBHs will eventually form a binary. 

Comparing the 1:2 mergers at the same Eddington 



fraction threshold yields insight into the AGN trigger- 
ing in the two simulations. At all Eddington fractions, 
the SMBHs spend more time active in the spiral-spiral 
merger. Luminosity thresholds favor detect ability in 
the elliptical-spiral merger because of the larger SMBH 
masses. We note that we have not included gas in the 
elliptical galaxy in the form of a hot halo. Ram pres- 
sure from the hot halo could strip gas from the spi- 
ral galaxy during the merger, preventing the secondary 
SMBH from accreting. Additionally, the hot halo could 
prevent stripped gas from cooling and reaching the pri- 
mary SMBH. We expect the AGN timescales from our 
elliptical-spiral merger, therefore, to be upper limits. 

1:10 Spiral- Spiral merger. We perform the same 
analysis on the 1:10 spiral-spiral merger with i g = 0.3 
from Callegari et al. (2011, cf. analogs of Fig. 1 for this 
simulation). This merger is longer than the 1:2 mergers 
and includes many more pericenter passages and stronger 
stripping of the secondary galaxy. The primary galaxy 
is relatively undisturbed, leaving the primary SMBH to 
grow quiescently. The secondary SMBH, on the other 
hand, accretes at a higher Eddington fraction until ram 
pressure strips all gas from the secondary galaxy. When 
the secondary galaxy enters the disk of the primary and 
its orbit circularizes, the secondary SMBH begins accret- 
ing again, leading to a phase of dual AGN activity at 
sub-kpc separations. Table 1 shows that the majority of 
dual AGN activity at high activity thresholds occurs at 
separations under 1 kpc. Overall, the 1:10 spiral-spiral 
merger yields less dual AGN activity than the 1:2 merg- 
ers, even though it is a factor of two longer in duration. 
Stronger tidal effects and ram pressure on the secondary 
galaxy mean that the secondary SMBH is unable to ac- 
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Figure 3. Evolution of the 1:2 elliptical-spiral merger. Panels 
and symbols as in Fig. 1. 

crete efficiently for a significant amount of time above 
kpc-scale separations. The AGN reach bolometric lumi- 
nosities of 10 43 erg s _1 for ~ 60 Myr and 40 Myr 
for the primary and secondary SMBHs respectively, and 
for ^ 7 Myr simultaneously. 

4. DISCUSSION AND CONCLUSIONS 

We perform numerical simulations of galaxy mergers, 
focusing on the separations and timescales for dual AGN 
activity. In all of the simulations presented here, much 
of the AGN activity is non-simultaneous at high activity 
thresholds. We find that the dual activity is generally 
a small fraction of the total AGN activity, with each 
SMBH accreting longer alone than in a pair. At our 
lowest thresholds, the long dual AGN timescales simply 
reflect the almost non-stop low level accretion onto the 
SMBHs. At higher thresholds we are instead probing 
accretion that is triggered by the dynamics of the merger, 
yielding better correlated accretion between the SMBHs. 
However, significant non-simultaneous accretion remains. 

Our results can be summarized as follows: 

• At high luminosity thresholds, almost all dual AGN 
activity occurs at separations < 10 kpc, where 
AGN triggering is strongest. Lower thresholds fa- 
vor instead large separations, where the galaxies 
spend more time during the merger. 

• Much of the AGN activity during the mergers is 
non-simultaneous. SMBHs are active alone rather 
than as a pair for ~ 90% of the time at L^i > 10 44 
erg s _1 (corresponding to an optical luminosity of 
order 10 43 erg s _1 , Marconi et al. 2004) , separa- 
tions > 1 — 10 kpc and Av > 150 km s _1 . These 
SMBHs will appear as either 'normal' single AGN 
(at low spectral or spatial resolution) or as offset 
AGN. This is in agreement with Comerford et al. 



(2009b), who find that it is more common to find 
offset AGN rather than dual AGN. 

• From Table 1, the expected dual AGN fraction with 
Lboi > 10 44 erg s _1 is 4-5% in our 1:2 mergers if 
we consider sufficiently luminous dual AGN that 
can be spectroscopically identified (Av > 150 km 
s _1 ). While our results are not directly compara- 
ble to observations of dual AGN, which typically 
occur at z < 1, the dual AGN fraction is in broad 
agreement with that found by Shen et al. (2011) in 
SDSS AGN. Less than 1% of luminous dual sys- 
tems would be identified via imaging in surveys 
with spatial resolution less than 10 kpc (Foreman 
et al. 2009). 

• The mass ratio and morphologies of merging galax- 
ies determine when their SMBHs will be fed and 
starved. Tidal stripping unbinds the stars and 
gas surrounding the secondary SMBH in our 1:2 
elliptical-spiral merger, preventing strong dual ac- 
tivity at sub-kpc separations. In the 1:10 spiral- 
spiral merger, the secondary galaxy is stripped of 
its gas during the fourth pericenter passage. 

We note that these estimate of dual AGN activity are 
upper limits as we have not included dust and obscura- 
tion in our models, nor dilution by star formation. This 
may hinder detection of the AGN (Sanders et al. 1988; 
Hopkins et al. 2005), especially in optical surveys, but 
perhaps not in hard X-rays (Koss et al. 2011). As we do 
not follow the SMBHs past our resolution limit, we do 
not include binary AGN or post-coalescence activity in 
our analysis. Additionally, there is likely to be consider- 
able variability in the SMBH acccretion rates on smaller 
timescales and spatial scales than we resolve (Hopkins & 
Quataert 2010; Levine et al. 2010). This unresolved vari- 
ability could decrease the overall dual fraction as more 
accretion becomes non-simultaneous. We will address de- 
tectability in different bands in a follow-up paper, taking 
into account the luminosity of the host and dust obscu- 
ration. 

The authors are grateful to the anonymous referee for 
comments that improved the clarity of the Letter. Simu- 
lations were run using computer resources and technical 
support from NAS. MV acknowledges support from SAO 
Award TM1-12007X, NASA award ATP NNX10AC84G 
and NSF award AST 1107675. 
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